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Abstract. We study a 2— parametric family of probability measures on the space of 
countable point configurations on the punctured real line (the points of the random 
configuration are concentrated near zero). These measures (or, equivalently, point 



processes) have been introduced in Part II (A. Borodin, math/9804087 ) in connection 
with the problem of harmonic analysis on the infinite symmetric group. The main 
result of the present paper is a determinantal formula for the correlation functions. 

The formula involves a kernel called the matrix Whittaker kernel. Each of its two 
diagonal blocks governs the projection of the process on one of the two half-lines; 
the corresponding kernel on the half-line was studied in Part III (A. Borodin and 



G. Olshanski, math/9804088) 



While the diagonal blocks of the matrix Whitaker kernel are symmetric, the whole 
kernel turns out to be J-symmetric, i.e., symmetric with respect to a natural indefi- 
nite inner product. 

We also discuss a rather surprising connection of our processes with the recent 



work by B. Eynard and M. L. Mehta (^ond-mat/9710230) on correlations of eigen- 
values of coupled random matrices. 



Introduction 

In this paper we continue the study of a remarkable family of stochastic point 
processes started in [01] (Part I), [Bl] (Part II), and [B01] (Part III). This Part 
IV is followed by the next paper [02] (Part V) which concludes our work on the 
subject. The summary of results from Parts I-V is presented in [B02]. 

The point processes in question live on the punctured interval / = [—1, 1] \ {0}. 
They are parametrized by two parameters z and z' which satisfy one of the following 
two conditions: 

(*) z' = z, z G C\Z; 

(**) z and z' are real and for a certain m e Z, m < z, z' < m + 1. 



We denote the corresponding process by V zz >. The definition of these processes can 
be found in [01, §4] and [Bl, Introduction]. 

The origin of our study is the work [KOV] on generalized regular representations 
of the infinite symmetric group — the processes V zz i govern, in a certain sense, the 
decomposition of these representations into irreducibles. 
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Part I was devoted to the translation of the representation theoretic problem in 
the language of the point processes, and to the calculation of the density function 
of the processes. In Part II we computed the higher correlation functions of V zz i 
and introduced a modification of the processes (so-called 'lifting') which substan- 
tially simplified the formulas for the correlation functions. Part III contains several 
applications of the results of Parts I, II. 

In Part II we gave explicit integral representations for all higher correlation func- 
tions, see [Bl, Theorem 1.7.3, Theorem 2.2.1]. However, these integral representa- 
tions were rather complicated. After lifting, see [Bl, Chapter 3], multidimensional 
integrals involved in the integral representations were reduced to two-dimensional 
ones. Moreover, we have proved that in the domains where all variables have the 
same sign, the correlation functions are given by determinantal formulas with a 
certain kernel on M + called the Whittaker kernel [B, Theorem 3.3.4]. 

In this paper we show that after lifting the correlation functions are everywhere 
given by the determinantal formulas with a certain kernel defined on M* . It is 
convenient to write this kernel in the matrix form according to the splitting M* = 
R+ U R_ . We call the kernel thus obtained the matrix Whittaker kernel. One of its 
blocks coincides with the Whittaker kernel introduced in [Bl]. 

The expression of the 'lifted' correlation functions via the matrix Whittaker 
kernel is our main result. 

Determinantal form for the correlation functions appears in many problems of 
random matrix theory and mathematical physics, see, e.g., [Dy], [Me2], [KBI]. In 
most situations the kernels are symmetric or hermitian (see, however, [B2]). But the 
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Remark 2.9 below. The appearance of kernels with such symmetry in determinantal 
formulas for correlation functions seems to be new. 

At the end of the paper we consider the systems of eigenvalues of two random 
coupled matrices studied in [IZ], [Mel], [MS], [EM], [MN], [Ey]. As was recently 
proved in [EM], the correlation functions of such systems are also given by deter- 
minantal formulas. We show that this result and our considerations have common 
combinatorial background. Here we give an alternative proof of the main result of 
[EM] in case of two coupled matrices ([EM] deals with a more general situation of 
several matrices coupled in a chain). Our proof seems to be more straightforward 
and less elegant than that from [EM]. 

The paper is organized as fol lows. In Section 1 we prove that the correla- 
tion functions of V zz > after lifting are given by determinantal formulas and express 
the kernel through two-dimensional distributions introduced in [Bl]. The proof is 
purely combinatorial. In Section 2 we show that the kernel obtained in Section 1 
can be conveniently written via the Whittaker functions. Section 3 deals with the 
eigenvalue correlations of two random coupled matrices. 

I am very greatful to G. I. Olshanski for his permanent attention and support. 
His 'jewel' interpretation of the proof of Proposition 1.1 made the presentation 
much more transparent. 



matrix Whittaker kernel turns out to be ' J-symmetric' where J 



see 



1. Determinantal formula 



Let us introduce the notation used in [Bl]. 
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Denote by § n ,d the set of mappings 



¥ >:{l,...,n}-{l,l , ;...;d,d'} 

subject to the two conditions 

1) ip is injective, i.e. (p(i) 7^ ip(j) if % 7^ j; 

2) Imy? fl {m, to'} 7^ for all to = 1, . . . , d. 

It is clear that Q n ,d 7^ if and only if n/2 < d < n. 

Starting from a function (or a distribution) F(r\, s\; . . . ; r^, Sd) in 2d variables 
we define the function (<pF)(xi, ■ ■ ■ , x n ) in n variables, cp G $ n ,d, as follows. Let us 
rename the variable by Xk if <f(k) = i, Sj by Xk if <fi(k) = j', and let us do this 
for all rj, Sj such that G Im</?. Then let us integrate F over all ri, s m such that 
/ ^ Im</?, to' ^ limp. The result is denoted by (<pF)(xi, . . . ,x n ). 

Let N(x,y) and w(x,y) be (generalized) functions in two variables such that 

supp N(x,y) C {(x,y) \x>0,y< 0}, 
suppty(x,y) C {(x,y) \x < 0,y > 0}. 

For all d = 1, 2, . . . set 



1 d 

# d (n,si;... ;r d , s d ) = ^[[^(s^r,) • det[iV(ri, Sj)]* j=1 . (1.1) 



i=l 



Our main result in this section is the following statement. 
Proposition 1.1. Let x\, . . . , x n 7^ 0. T/ien 

(p#d)(si, •••,**) = det[iV(x i ,x J -)]"j=i 



(1.2) 



where 



N(x,y) 



J N(x, s)w(s, y)ds, 



x, y > 



iV(x, y), x > 0, y < 

JJ w(x,r)N(r, s)w(s,y)drds — w(x,y), x < 0, y>0 



J" w(x, r)N(r, y)dr 



x, y < 



Remark 1.2. Theorem 3.2.6 of [Bl] claims that the nth lifted correlation function 
of our processes is equal to the left-hand side of (1.2) for the following choice of 
N(x,y) and w(x,y): 

N(x, -y) = t ( ( p z/ (x) ( P- z/ (y)e- x - y )Q((p z ,- z -i(x) ( P- z , z - 1 (y))Q(S(x - y)x[o,i](x)) , 

(1.3) 

w(x,y) = TT —, (1.4) 
\x\ +y 

see Section 1.4 and 1.6 of [Bl] for the notation. 
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(Note that our function N here differs from N defined by [Bl, (3.6)] by changing 
the sign of y and multiplication by the constant t = zz' '.) 

Proof of Proposition 1.1. The proof will consist of five steps. 

Step 1. On the set <& n ,d there is a natural action of the symmetric group Sd 
which permutes the pairs (i, i') for i = 1, . . . , d. On the other hand, the function 
Hd defined in (1.1) is invariant with respect to the permutations of the pairs (r i? s^), 
i = 1, . . . , d. This implies that all elements <p G <& n ,d lying in the same orbit of 
Sd give the same contribution to the left-hand side of (1.2). Each orbit contains 
dl elements, hence, we may assume that the summation in the left-hand side of 
(1.2) is taken over the set of orbits of Sd in <& n ,d and throw out the factor in the 
right-hand side of (1.1), cf. [Bl, Proposition 1.7.1]. 

Step 2. Let us plug the formula (1.1) into the left-hand side of (1.2) and expand 
all the determinants. Then the summands will be parametrized by the triples 
(d, <p, a) where <p is an orbit of Sd in $ n> d and a G Sd parametrizes the expansion 
of the determinant (1.1). We may also consider a as a map from {1,2,... , d} to 
the set {l',2', ... ,d'}. 

We associate to each such summand an oriented (not necessarily connected) 
necklace with 2d black and white beads, some of the beads are marked by the 
numbers from 1 to n, in the following way. 

Let us fix a representative ip G (p C $ n ,d- 

The necklace contains d white and d black beads, every two neighbors are of dif- 
ferent colors. White beads correspond to the elements {1,2,... , d} and black ones 
- to the elements {1', 2', . • • , d'} of the set {1, 1'; . . . ; d, d'} used in the definition of 

A white bead % is marked by the number k if (p(k) = i; a black bead j' is marked 
by the number / if ip(l) = j' . 

The white bead i always follows the black bead i' and is followed by the black 
bead o~(i). 

Now we forget that the beads correspond to the elements of {1, 1'; . . . ; d, d'} and 
obtain an oriented necklace with black and white beads some of which are marked 
by 1,2, ... ,n. 

It is not difficult to see that different representatives of the same orbit <p produce 
the same necklace, and different orbits correspond to different necklaces. 

For a fixed necklace the corresponding summand of the left-hand side of (1.2) 
can be constructed as follows. We take 2d indeterminates and associate them to 
the beads somehow. Then for each edge (interval between two beads) we produce 
a factor like this: 

for an edge going from a white bead (indeterminate r) to a black bead (indeter- 
minate s) we take N(r, s); 

for an edge going from a black bead (indeterminate s) to a black bead (indeter- 
minate r) we take w(s, r). 

Thus, we obtain d factors, we multiply them and attribute to the indeterminates 
corresponding to the beads marked by 1, . . . , n the names x±, . . . ,x n respectively. 

Finally, we integrate over all other (mute) indeterminates and add the sign sgna. 
The sign will be discussed in Step 5. 

Step 3. Now we shall translate to the 'jewel language' the right-hand side of 
(1.2). 

We may assume that all variables xi, . . . ,x n have definite signs. 
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First, we expand the determinant in the right-hand side of (1.2), the summands 
are parametrized by permutations r G S n . Further, if a summand contains a factor 
N(x, y) with negative x and positive y, then 

N{x,y) = J J w(x,r)N(r, s)w(s,y)drds — w(x,y), (1-5) 

r,s 

and we break this summand into two: one will contain 

//»(*,rWr, .)»(.,»)**, (1.6) 

r,s 

and the other 

-w(x,y). (1.7) 

For each summand we produce an oriented (not necessarily connected) necklace 
with n beads corresponding to the variables xi, . . . ,x n ; the beads corresponding 
to positive variables are white, to negative - black. All beads are marked by the 
numbers 1, . . . , n according to the numbers of corresponding variables. The bead 
marked by i is followed by the bead r(i). 

On each edge going from a black ('negative') bead i to a white ('positive') bead 
j we put the sign W or according to which of the summands (1.6) or (1.7) was 
chosen while breaking N(xi,Xj) into two parts. 

In the first case we shall say that we used iV-choice, in the second - w- choice. 

Step 4. Now we shall establish one-to-one corresponds between the sets of 
necklaces introduced in Steps 2 and 3. 

First, let us take a necklace from Step 2. It contains 2d black and white beads 
some of which are marked by the numbers 1,... , n. Then we can immediately 
produce an oriented necklace with n beads by throwing out all mute (unmarked) 
beads. 

The second condition in the definition of Q n ,d implies that at least one end of 
each edge going from black to white is not mute (at least one element in each pair 
is covered by ip G $ n ,d)- It means that after removing all mute beads, we 
can obtain an edge 'black-to-white' in two different ways: it either comes from the 
initial necklace unchanged, or there were two mute beads removed from the middle 
of this edge. In the first case we put the sign l w : on this edge, in the second - the 
sign 'AT. 

Thus, starting with a necklace from Step 2 we produce a necklace from Step 3. 
Let us go the other way around. 

We take a necklace with n beads from Step 3. Then we follow the rules: 

(1) inside any edge of the type 'white-to-white' we insert one mute black bead, 

(2) inside any edge of the type 'black-to-black' we insert one mute white bead, 

(3) the edges 'white-to-black' remain unchanged, 

(4) the edges 'black-to-white' with the sign remain unchanged, 

(5) inside any edge 'black-to-white' with the sign 'AT we insert one white and 
one black mute beads, white goes first. 

Clearly, two procedures described above establish one-to-one correspondence 
between the necklaces of Steps 2 and 3, thus, between the summands of the left- 
hand side and the right-hand side of (1.2). 
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Straightforward check shows that up to sign the corresponding summands from 
Steps 2 and 3 are equal. To compute, up to sign, the summand corresponding to 
a pair of necklaces, we take each edge of the necklace with n beads and according 
to 5 different situations for this edge described above (suppose the edge goes from 
the zth bead to the jth bead) we produce the factor: 



(1) : / N(xi, s)w(s,Xj)ds; 

J s 

(2) : / w(xi,r)N(r, Xj)dr; 

J r 



(3) : N(xi,Xj); 

(4) : w(xj,Xi); 



(5): 



w(xi, r)N(r, s)w(s, Xj)drds — w(xi, xj). 



r,s 



Then all such n factors are multiplied. 

Step 5. Let us verify that corresponding summands have the same sign. 

We take two corresponding necklaces: one from Step 2, the other from Step 3. 

The sign for the first one equals sgna (a E Sd was introduced in Step 2). The 
number of cycles of a is equal to the number of connected components of the 
necklace, let us denote it by c. Then sgncr = ( — l) d ~ c . 

The sign for the second necklace is equal to sgnr (r was introduced in Step 
3) times —1 to the power of the number of w- choices, because of the sign ' — ' in 
(1.7). The number of cycles in r is equal to c, because both necklaces have the 
same number of connected components. Thus, sgnr = (— l) n_c . Furthermore, the 
number of w-choices is equal to n — d, it can be easily verified by induction on the 
number of w-choices. Hence 

n — c — number of w-choices = d — c, 

and the signs are equal. □ 

2. The matrix Whittaker kernel 

Our goal in this section is to compute explicitly the kernel N(x,y) introduced 
in Proposition 1.1 for the special N(x,y) and w(x,y) given by (1.3) and (1.4). 
It is convenient to write N(x,y) in the matrix form 



N-+(x,y) N—(x,y) 



where 



N++(x,y) 

N + _(x,-y) 
N_+(-x,y) 



J N(x,s)w(s,y)ds, 

s 

N(x,y), 

w(x, r)N(r, s)w(s, y)drds — w(x, y), 



r,s 



N {—x,—y)= / w{x,r)N{r,y)dr, 



x,y>0; 

x > 0,y < 0; 
x < 0,y > 0; 

x, y < 0. 



Starting from now we set, see Remark 1.2, 

N(x, -y) = t (cj) z/ (x)(j)- z/ (y)e- x - y ) Q(</> z ,-z-i(x)<l>- z , z -i(y))® (6(x - y)x[o,i](x)) , 

(2.1) 

f n^— , x < 0, y > 
= <^ lxl+y . (2.2) 

t 0, otherwise 

Moreover, in all our calculations below we shall assume that 

-1 < 9fcz, Sfcz' < 0. (2.3) 

Under this assumption several distributions used below become just integrable func- 
tions, which will be suitable for us. Our final result will depend on z and z' ana- 
lytically. Thus, by analytic continuation, it will hold for all pairs z, z' . 

We use the standard notation W Kjfl (x) for the Whittaker functions, see [E, Chap- 
ter 6] . Note that 

W^{x) = W K ,^(x). 

Proposition 2.1. 

z + z' 

, , , T . , sin Ttz sin tyz' ( x\ 2 -*- v , , i 
N+-(x,y) = N(x,-y) = = ~ e 2 M 2 

W z+Z ' +1 ,z-z> (x)W -z-z'+i ,z-z> (y) + tW z+z'-i ,z-z> (x)W - z - z '-i , z -z' (y) 

^ 2 '^2 2 '^2 2 2 2 2 

x + y 

(2.4) 

Proof. Let us multiply (2.4) by (x+y) and compare the (fc, /)-moment of both sides. 
We shall need the relations 

(M*)e- X , x m ) = (a + l) m , (0„, 6 (x), x m > = r(fl 2) 

which are just Euler gamma and beta integrals; and 



/ 



2 a '2 2 V ' r(a + 6-c + l) 



which is a version of [E, 6.10(7)]. 

The moments of a pseudoconvolution of distributions are the products of mo- 
ments of the factors. Thus, the (k, /)-moment of (x + y)N(x, y) equals 

((x + y)N(x, -y), *V> = t ^±A^' + + + 1)« 



+ * 
= t 
+ t 



(k + l)W(k + I + 2) 
(z + l) k (z' + l) k (-z + l)i+i(-z' + 

k\{l + l)\{k + I + 2) 
(z + l) k (z' + l) k (-z + l)i(-z' + 

k\l\ 

2 (Z + I) k (z' + I) k (-Z + I) l (-Z' + 1); 

(fc+l)!(Z + l)! 
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Straightforward computation shows that these two summands are exactly the 
(k, I) -moments of two summands in the numerator of the right-hand side of (2.4). 
Thus, we proved that the moments of both sides of (2.4) coincide after multiplication 
by (x + y). 

In general, it is not sufficient to conclude that our distributions coincide, because 
they are not compactly supported. However, it is easy to check that both sides of 
(2.4) after multiplication by (x + y) can be realized as liftings of some compactly 
supported distributions. (See Chapter 3 of [Bl] for the definition of the lifting.) By 
comparison of moments, these compactly supported distributions must coincide, 
and so do their liftings. □ 

Proposition 2.2. 



/ \ f N(x,—s) , 1 { x\ 2 y-x , .1 

N++(x,y) = / ; U s = - e—(xy) * 

J s + y T{z)T(z') \yj 

s 

W z+Z > +1 , z - z > (x)W z+z >-! , z - z > (y) - W z+z '-i , z - z > (x)W z+z ' +1 , z - z > (y) 

^ 2 '_ 2 2 '^2 2 '^2 2 >_ 2 

x-y 



(2.5) 



Proof. This claim is implicitly proved in [Bl]. Let us explain why. 

Application of the Stieltjes transform to a pseudoconvolution is equivalent to the 
application of it to one of the factors, if the others does not have a singularity at 
zero. Assuming that —1 < $tz' < 0, see (2.3) above, we may take the Stieltjes 
transform of the second factor of (2.1), instead of applying it to the whole N(x, y). 

Then we use Lemma 2.2.3 from [Bl] which claims that 



y- z (l + yy-'. (2.6) 

s + y 

s 

We get 

N ++ (x,y)= [ N{x '~ s) ds = t(<p z ,(x) ( j ) -Ay)e- X - V ) 
J s + y 

s 

y~ z (l + yY' 1 ) (S(x - y) X [o,i](x)) . 

By definition of pseudoconvolution (see [Bl, 1.4]), we get 

N ++ (x,y) = t J <j> z ,{u)<l>_ z ,{v)e- u - v 

u,v 

(2-7) 

// x \ ( y \~ z ( ' y \ z ~ 1 dwdudv 
4>z,-z-l — — ! + — — . 

\uw / \vw/ V vw / w z uv 



X 





cf. [Bl, (2.6)]. 

Lemma 2.2.4 of [Bl] states that 



J \uwJ \vw / \ vw/ 



r 2 \ z ~ 1 dw 4>- z (u - ri)(/> 2 (r 2 + v) 



uvw 2 r\v + r-iu 

o 



Applying this to (2.7) we get 

/ e -u-v 
(j) z/ (u)(j)- zl (v)(j)- z (u - x)4> z {v + y) ■ dudv. 
xu + yv 



By introducing new integration variables 

ti = u/x — 1, t 2 = v/y 

we see that 



z 

X \ y—x 



N++(x,y) = tl-J e~ 

ff e -x(t 1 + l/2)-y(t 2 + l/2) 

x 11 0_ a (ti)0_^(t 2 )^'(ti + l)<f>z{ti + 1) 1 +t 2 + l dtldh ' 

tl ,*2 

cf. [Bl, Theorem 3.3.1]. The formulas (3.8) and (3.9) from [Bl] (together with the 
definition of the kernel M(x, y), see [Bl, Theorem 3.3.1]) conclude the proof. □ 

By changing the signs of the parameters z and z' in Proposition 2.2 and changing 
x <-> y, we get the following statement. 

Corollary 2.3. 

N—{x,y) = [ N{r '~ y) dr = 1 (-) 2 e^ixy)^ 

J r + x r{-z)T(-z') \yj 

r 

W - z - z > +1 +z - z > (x)W - z - z '-i , z - z > (y) - W -z-£-i +z - z > (x)W - z - z >+i +z - z > (y) 

^ 2 '^2 2 '^2 2 2 2 2 



x-y 

(2.8) 



Finally, we compute N y(x,y). 

Proposition 2.4. 



N(r,-s) 1 



iV_+(a;,y) = / / -———-drds — = - - e » (xy) * 

(r + x)(s + y) x + y \y/ 

± *-s (x)W z+z>+1 , (y)+tW - z - z ,.i LZ .Ax)W z+z/ . 1 ± w (y) 

^ 2 ' 2 2 ' 2 2 ' 2 2 ' 2 

x + y 

(2.9) 

Remark 2.5. Note that 

N_ + (x, y) = —. - e x+ VN+-(y, x). (2.10) 

SID. It Z SIDTIZ' 



This is a non-trivial and surprising fact. A detailed discussion of this 'coincidence' 
(or rather of its modification, see Remark 2.9) can be found in [02]. 
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Proof of Proposition 2.4- First, as in the proof of Proposition 2.2, we apply both 
Stieltjes transforms to the second factor of the pseudoconvolution (2.1). Using (2.6) 
twice, we get 

/ *( r >-*l drds = t^A^-Av)e- x - y ) 
J {r + x){s + y) v 

r,s 

(x z (l + x)-*- 1 y-*(l + y) z - 1 ) (S{x - y)x[o,i]{x)) . 
By definition of the pseudoconvolution, we get 

N ( r >~ s > -drds = t J (P z ,(u)(P- z/ (v)e- u - v 

u,v 

1 

f ( X V (l X X~ Z ~ X ( y \~ z (i y \ z ~ 1 dwdudv 
J \uwJ V uwJ \vwJ V vwJ 



(r + x)(s + y) 

u,v 

1 



X 



W 2 UV 







Introducing new variables 



we obtain 



T\ = u/x 
r 2 = v/y 



Tl,T 2 

1 



x J (1 + r 1 w)~ z ~ 1 (l + T 2 w) z ~ 1 dwdT 1 dT 2 . 



We shall need the following lemma. 
Lemma 2.6. 



l 

/ P + ^>->-Hi + w) -<to = ((1±2) - 1) 



Proof of Lemma 2.6. First we use [E, 5.8.2(5)]: 
l 

J(l + T lW )- z -\l + r 2 w) z - 1 dw = Fi(l, ^ + 1,-^ + 1, 2; -n - r 2 ), 
o 

where Ft (a, /?', 7; x, y) is the Appell's hypergeometric function. By [E, 5.10(1)], 
we can simplify the last expression: 

Fi(1,z + 1,-z + 1,2;-ti-t 2 ) = (1 + t 2 )- 1 F ^ z + l;2;^-^j . 
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But 

F(M+l;2;r) = — ((1-t)"*-1), 



ZT 



and we arrive at our claim. □ 
Applying the lemma, 



Tl,T 2 



(2.11) 



(.r + //)< " " 7 - ! ( tt" - " ~rr~ ) e" xri 



Let us multiply this relation by x + y. Note that 

d_ 

dr\ <9t 2 

cf. proof of [Bl, Theorem 3.3.4]. 

We are going to integrate (2.11) multiplied by x + y by parts. The following 
formula is easy to check 

d d \ f 1 1 / / 1 + r 2 x 



+ 



dr\ <9r 2 



Z T 2 -T! V V 1 + Tl 



11 

- <f>z'-l(Tl)<f>-z'-l{T2) ~ - z '-i(ti)0_^_i(t 2 )(1 + Ti) Z (1 + T 2 )* 



0^(r 1 )0_^(r 2 )(l + t x )-*-\\ + r 2 )-- - 1 



Hence, integrating by parts, 



+ / ? ^T7 ~ , drds= (-) I ( z /_i(ri)0_ z /_i(r 2 ) 



(r + x)(s + y) \j/ 



Tl,T 2 



- ^_ 1 (r 1 )^_ 1 (r 2 )(l + n)-*(l + r 2 ) 2 
-t^(r 1 ) ( />_ 2 ,(r 2 )(l + r 1 )- 2 - 1 (l + r 2 ) 2 - 1 ^e-^-^rfr 1 rfr 2 . 

Using Euler gamma integral for the first summand and the standard integral rep- 
resentation of the Whittaker functions 

W K ^(x) = e-^x^ 1 ' 2 J M _ K _ 1/2 (r)(l + rr+^e-^dr 

T 

for the last two, we rewrite the last relation in the form 

z + z' 

N(r, —s) fx* 



(x + y) / r- rdrds = 1 — — e 2 (xy) 2 

! (r + x)(s + y) \y> 



r,s 



x ( W - z - z '+i ,' -z {x)W z+z'+i £ _^(y) + tW -z-z'-i ,' -z (x)W z+z>-i , z > (y) 

' 2 ' 2 2 ' 2 2 ' 2 2 ' 2 
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which is equivalent to (2.9). □ 

Note that in all four quadrants the expression for N(x, y) has the same factor 



e 2 



\x\ 2 e 2 

i i z + z ' _ v 

\y\ 2 e 2 



(2.12) 



(For example, if x, y < 0, N(— x, —y) = iV__(;r, y) contains the factor 



x 



y 



z + z' 
2 



which after changing the signs of x and y coincides (2.12).) 

But all factors (2.12) disappear in the determinants of the form det N(xi, Xj). 
Thus, we can introduce a new kernel as follows. 

Denote 

A+(x) = x~^W z+z>+i + z-z< (x), 

2 '-"-2 

B+(x) = x~^W z+z*-i j_z-z> (x), 

2 > 31 2 



A_(x) = x 2 W -z-z'+i 



B_(x) =x-*W- 



z — z 1 — 1 I z — z 

2 2 



.(*), 
.(*)• 



For x, y > set 

K(x,y) = K ++ (x,y) 
K(x,-y) = K+_(x,y) 
K(-x,y) = K_+(x,y) 



A + (x)B + (y) - B + (x)A + (y) 



x - J/ 



Vsin7r2 sinyr^' A + (x)A_(y) + t B + (x)B_(y) 
n x + y 

Vsin 7T^ sin rcz' A + (y)A_(x) + t B + (y)B_(x) 



K{-x,-y) = K__{x,y) = 



7T 
1 



T(-z)T(-z') 



x + y 

A_ (x)B_(y) - B_(x)A_(y) 
x-y 



or, in matrix form, 



K{x,y) 



K++(x,y) K+-(x,y) 



We call K(x,y) the matrix Whittaker kernel 

All our previous work is summarized by the following statement. 

Theorem 2.7. The nth correlation function of the lifted process V zz > has the form 

Pn ^ { x li • • • > x n) = det[K(xi, Xj)]™j—i 

where K(x,y) is the matrix Whittaker kernel. 
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Remark 2.8. Theorem 3.3.4 of [Bl] where we established the determinantal form 
of the lifted correlation functions on a part of the phase spase is a direct corollary 
of Theorem 2.7. The Whittaker kernel introduced in [Bl] is the '++' block of the 
matrix Whittaker kernel. 

Remark 2.9. Note that all blocks of the matrix Whittaker kernel are real-valued 
functions and 

K++(x,y) = K++(y,x), K—(x,y) = K—(y,x), 
K+-(x,y) = -K_+(y,x), 



cf. Remark 2.5. This means that the matrix Whittaker kernel is not symmetric 
but ' ./-symmetric' for 

"1 ' 

-1 



J = 



A general discussion of stochastic point processes governed by such kernels can be 
found in [02]. 



3. TWO RANDOM COUPLED MATRICES 

In this section we will explain the connection of our problems with the system 
of eigenvalues of two random coupled matrices. This system was introduced in 
[IZ] and studied in [Mel], [MS], [EM], [MN], [Ey]. In particular, B. Eynard and 
M. L. Mehta proved in [EM] that the correlation functions of this system are given 
by determinantal formulas (in fact, they proved the determinantal formulas for the 
correlation functions in case of finitely many matrices coupled in a chain). Let 
us describe this result in more details. We will follow [EM] slightly changing the 
notation. 

Consider two complex hermitian N x N matrices A and B with the probability 
density 

F(A, B) = const ■ exp(- tr{U(A) + V(B) + cAB}) 

where U(x) and V{x) are real polynomials of even degree with positive coefficients 
of their highest powers and c is a real constant. 

Let us denote the sets of (real) eigenvalues of A and B by {xn, . . . , £iat} and 
{a?2i, • • • , %2n}- As was proved in [IZ], [Mel], the probability density for the eigen- 
values has the form 

Pn(xii, • • • , x 1N ;x 2 i, • • • , x 2 n) 

= const- det[w(x 2 i,x lj )]f I j=1 ■ JJ[(x H - x lj )(x 2i - x 2j )] ( 3-1 ) 

i<j 

where 

w(x, y) = ex.p(-U(y) - V(x) + cxy). 
The correlation functions are defined as follows: 

Pk,i{xn, . . . ,x lk ; x 21 , ... , X21) 

AH 2 f , , , (3.2) 



(N-k)l(N-l) 



y J Pn{xi, x 2 ) dx lk+1 . . . dx 1N dx 2i+1 . . . dx 2N . 
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Under certain non-degeneracy conditions on the pairing 

(f(x),g(y)) = J f(x)g(y)w(y,x)dxdy, (3.3) 

we can find two systems of biorthogonal polynomials 

{P (x),... .PN-tix)} and {Qo(y),... ,Q N -i(y)} 
with respect to this pairing. In other words, deg Pi = deg Qi = i and 

{Pi i Q j ) $ij ■ 

Set 

N-l 

H(x,y)= J2 P i( x )Qi(v) ( 3 - 4 ) 

i=0 

Proposition 3.1 ([EM]). 

Pk,l(xil, ■ ■ ■ ,Xlk',X21, ■ ■ ■ ,X2l) = det[Kij(x ir ,Xj s )]ij =li 2;r=l,...,k;s=l,...,l (3.5) 

where 

Ku(x,y) = J H(x,s)w(s,y)ds; 

s 

K 12 (x,y) = H(x,y); 

K2i(x,y) = J J w(x,r)H(r, s)w(s,y)drds — w(x,y); (3-6) 



K 2 2(x,y) = J w(x,r)H(r,y)d7 



Remark 3.2. Note that the formula (3.6) is identical to the formula for the kernel 
N(x,y) in Proposition 1.1. This is not an accidental coincidence. It turns out that 
there exists a wide class of measures on the infinite-dimensional Thoma simplex 
such that after lifting the correlation functions of these measures are given by de- 
terminantal formulas with a matrix kernel. Moreover, all these matrix kernels have 
the form (3.6) for an appropriate distribution H(x,y) (e.g., the matrix Whittaker 
kernel) . 

In particular, this class of measures include many finite-dimensional measures 
for which the techniques of [EM] can be applied. For example, one of such measures 
after lifting has the density of the form (3.1) with the weight function w(x, y) = -^^j 
on R + x R + , cf. (2.2). 

Now we shall give a sketch of the proof of Proposition 3.1 which is different from 
the proof given in [EM]. Namely, we shall reduce the statement to Proposition 1.1 
proved in Section 1. 

Sketch of the proof of Proposition 3.1. Following the notation of Section 1, let us 
introduce the set &k,i;d of mappings 

(p : {il, . . . , i/c; 2 1, • • • , 2^} -> {1, 1'; • • • , d, d'} 
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subject to the three conditions 

1) <p is injective; 

2) Imp fl {m, m'} 7^ for all m = 1, . . . , <i; 

3) (/j(im) G {1, . . . , d}, p(2Tn) G {1', • • • , d'} for all m. 

By analogy with Section 1, for each <p G and a function F(ri, si; . . . ; rd, Sd) 

in 2<i variables we define the function (<pF)(xu, . . . , x±k', X21, • • • , a^z) in /c + Z vari- 
ables. Specifically, according to ip we attribute to some variables ri names xi m , to 
some variable Sj the names x^m, and then integrate F over all mute variables. 

We shall prove that 

Pk,l(xn,... ,X!k-,X21, ■ ■ ■ ,X2l) = ^ ^ ((pHcfiixu, . . . ,X 1 k-,X21, ■ ■ ■ ,X2l) 

d <p€&k,l;d 

(3.7) 

for 

1 d 

i7 d (ri,si;... ;r d ,s d ) = — JJ^Sj,^) • det[H(n, Sj)]f j=1 . (3.8) 

' i=l 

Then word-for-word repetition of the 'jewel' proof of Proposition 1.1 will prove 
(3.5); the only difference is that in this case we consider the necklaces where white 
beads are marked by the elements of the set {jl,... , ik} and black beads are 
marked by the elements of the set {2I, . . . , 2^}- 

In order to prove (3.7), we shall first show that the sets ^k,i-d naturally appear 
in the expansion of the determinant of a matrix along k rows and / columns and 
then apply this observation to the determinant in the RHS of (3.1). 

Let M = (iriij)fj =1 be a sufficiently large square matrix. For any (p e &k,i-,d we 
define (<pM) as a certain sum of expressions of the form 

±m aibl ■■■m adbd M^ ^ y ( 3 - 9 ) 

where by M(''^'''? d we denote the determinant of the submatrix of M obtained 
from M by removing the rows ai, . . . ,ad and the columns b±, . . . ,bd- 

The sum is defined as follows. We associate the indices {a±, ... , ad', 61, . . . , bd} 
with the elements of {1, ... , d; 1', . . . , d'} and the numbers {1, . . . , k; 1, . . . , 1} with 
the elements of {il, . . . , ±k; 2I, • • • , 2^}- Then, according to our map 

ip : {il, . . . , ifc; 2 1, • • • , 2/} -> {1, • • • , d; 1', . . . , d'} 

we specify some of a^s and fej's. As for the rest, we let d—k other (not specified) a^s 
vary over all (d—k) -tuples of pairwise distinct numbers from the set {k+1, . . . , iV}, 
and d — I other b^s vary over all (d — /)-tuples of pairwise distinct numbers from 
{l + l,...,N}. 

The sign in (3.9) is chosen in such a way that (3.9) enters the expansion of det M 
with positive sign. 

Lemma 3.3. For any integer k, I > and sufficiently large square matrix M 

d ¥>€*k,j ; d 
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This formula is the expansion of det M along the first k rows and / columns. For 
example, for k = 1, I = 0, (3.10) turns into the usual first row expansion of detM: 



N f 1 \ 

detM = ^(-l)* +1 m H -Ml , J. 
i=i \---i---y 



For k = 1, / = 1 we get the formula 

det M = mil • M (I " + ^ (-l) l+i+1 m i imi J - M^j "j. 

Here the first term corresponds to the unique element of 3>i,i ; i and the second term 
corresponds to the two elements of $i,i ; 2, which give the same contribution. 

For arbitrary k and / = 0, (3.10) coincides with the well-known Laplace expan- 
sion of the determinant. 

The proof of Lemma 3.3 is straightforward. The factor in (3.10) may be 
removed if the summation is taken over the sets of orbits of S^'s in <&/<;, cf. Step 
1 of the proof of Proposition 1.1. 

Using elementary row and column transformations in the Vandermonde deter- 
minants of (3.1), we can rewrite (3.1) in the form 

PivOm, . . . ,x 1N ; x 2 u ■ ■ ■ ,x 2 n) 

1 (3.11) 

= ^• det H^2i,£ij)] i , J= i • &et[Pi{x l0 )] i73=1 ■ [Qi(x 2j )] ijj=1 . 



Now we apply (3.10) to the first determinant of (3.11). We shall show that after 
normalization and integration over extra variables as in (3.2), the </?-term of the 
expansion of (3.11) will give exactly the 99-term of the RHS of (3.7) for all <p's. 

So, we set M = (w(x2i, x ij))fj=\ an d employ Lemma 3.3. 

Note that for a fixed ip the contributions of all summands (3.9) to the RHS 
of (3.2) are equal because they differ by permutations of the integration variables 
Xik+i, ■ ■ ■ ,xin and X21+1, • • • , x%n- The total number of summands (3.9) is the 
number of {d — A;)-tuples and {d — Z)-tuples of pairwise distinct numbers taken from 
the sets with N — k and iV — / elements respectively. This number equals 

(N-k)l (N-l)l 

(N-d)\ (N-d)V [ ' 

So we keep this combinatorial factor and choose one suitable for us summand of 
the form (3.9), namely, such that 

{ai,... ,a d } = {l,...,d}, 
,b d } = {l,... ,d}. 

Then this summand has the form 

d 

^ 2 -Y[w(x 2t(i) ,x u ) ■ det[«;(x2i, a;ij)]5=d+i • det i p *( x ij)}Z=i ' [Qi( x 2j)}i*,j=i- 

(3.13) 
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for a certain r G S d - 

The expression (3.13) can be easily integrated over xid+i, ■ ■ ■ £;uvj x 2 d+i, ■ ■ ■ , %2N- 
By expanding all the determinants, using orthogonality conditions and Gram's 
formula we get 



(N — d)l 2 d 

y^— Ylw(x 2T (i),x u ) ■ det[H(x 1 i,x 2 T{i))}i t j= 1 



(3.14) 



where H (x, y) is defined in (3.4). Finally, integrating (3.14) over xik+i, ■ ■ ■ , X\d an d 
X21+1, ■ ■ -X2d and multiplying it by (3.12), combinatorial factor from (3.2) and 4 
from (3.10), we arrive at (ipH d )(xu, . . . , xi k ; x 2 i, ■ ■ ■ , x 2 i). □ 
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